Possible Realization of non-BCS type Superconductivity in gapped Graphene 
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We show that the gauge field induced due to non-uniform hopping, in gapped graphene, can 
give rise to a non-BCS type of superconductivity. Unlike the conventional mechanisms, this su- 
perconductivity phenomena does not require any pairing. We estimate the critical temperature for 
superconducting-to-normal transition via Berezinskii-Kosterlitz-Thouless mechanism. Possibility of 
observing the same in ultra cold atomic gases is also pointed out. 
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Graphene is an atom thick allotrope of carbon, which 
has unusual electronic properties, owing to its hexagonal 
honeycomb lattice structure. It was first shown by Wal- 
lace that, the valence band and conduction band touch 
each other at six points in the momentum space, around 
which quasiparticle dispersion is linear Later, Se- 
menofF demonstrated that, the low energy quasiparticle 
excitations around these points, out of which only two are 
inequivalent, actually satisfy the massless Dirac equation 
Here, the Lorentz symmetry is an emergent feature, 
due to the honeycomb structure of the material, rather 
than a true space-time symmetry, and the speed of light 
here is replaced by Fermi velocity of the electrons. Ex- 
perimental isolation of this material opened a way to sim- 
ulate relativistic quantum physics and check validity of 
various fertile ideas [3]. Infact, one of the very important 
consequences of relativistic quantum physics, Klein para- 
dox, was experimentally seen for the first time, certifying 
the validity of the underlying physics 0] . Amongst other 
interesting effects, that make this material important, are 
room temperature quantum Hall effect J5|, Andreev re- 
flection 0, and universal conductance 0, Q- Graphene, 
being intrinsically a two dimensional system, can lead to 
realization of peculiar effects, only possible in the pla- 
nar world, like fermion number fractionalization due to 
vortices [9| and anyonic excitations |1C| . 

Possibility of superconductivity in graphene has at- 
tracted considerable attention from both theoretical and 
experimental side [ill However, the search for su- 
perconductivity in pure graphene has yielded null re- 
sults, till this date. Proximity induced superconductiv- 
ity in graphene has been experimentally seen, indicat- 
ing that graphene has phase coherence [Tj]. The possi- 
bility of superconductivity in pure and doped graphene, 
and graphitic layers, has been theoretically investigated, 
using variants of conventational pairing mechanisms 
[3 [l^. In Ref. [16], the authors address the possi- 
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bility of high Tc superconductivity using resonanting va- 
lence bond model for doped graphene. Doped graphite is 
known to be superconducting, and possible mechanisms 
giving rise to this effect have been proposed [13, [3 ■ In 
Ref. [13 , the authors propose a mechanism of supercon- 
ductivity arising due the edge states in graphene. 

In this letter, we show that, apart from the conven- 
tional BCS type pairing mechanisms, two dimensional 
nature of graphene can make it superconducting via 
an unconventional non-BCS type mechanism. We show 
that, this mechanism does not require any short range 
four fermi interaction usually associated with pairing. 
Rather this mechanism crucially relies on presence of a 
gauge field, which may arise from lattice deformation, 
which locally modify the hopping amplitudes. It is shown 
that, existence of such a gauge field interacting with mas- 
sive Dirac fermions, where mass may arise due several 
plausible mechanisms, via one loop correction can cou- 
ple to external photon field through a topological Chern 
Simons coupling. This topological coupling induces a 
parity preserving non-topological mass for the photon 
field, giving rise to perfect diamagnetism or Meissner ef- 
fect, which is a hallmark of superconductivity. Further, 
planar nature of this model ensures logarithmic confine- 
ment of all the vortex excitations present in the induced 
gauge field, untill a critical temperature, above which the 
via Berezinskii-Kosterlitz-Thouless mechanism, the vor- 
tex unbinding takes place and superconducting state is 
destroyed. We estimate the critical temperature for this 
model and show that it is proportional to the gap. We 
conclude by pointing out a possible way of realizing the 
above phenomenon in ultra cold atomic gases, in optical 
lattice with tunable interactions. 

In tight binding approximation, the electron hopping 
Hamiltonian, defined on a hexagonal lattice with hopping 
energy ii(r), reads: 

H = -Y.Y1 U{v){o){v)b{v + s,) + b\v + s,)a{^)), 

r 1=1,2,3 

(1) 

where the fermion operators a and b act on the two in- 
terpenetrating sublattices A and B [20]. The hopping 
energy may be written as ii(r) = t + Sti{r), where t is 
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constant, so that the above Hamiltonian can be written 
conveniently as: 



— 6_, a_, a_(-)-^, so that the full Hamiltonian 



H = Hd+Hg 



(2) 



where, 



Hd^-Y, i(a^(r)6(r + Si) + 6t(r + Si)a(r)), 

r i=l,2,3 

Hg = -Y Y '5<.(r)(at(r)6(r + Si) + 6t(r + Si)a(r)) 

r 1=1,2,3 

As is well known, that can be linearized around two 
Dirac points K±, to yield: 



Hr 



J (i2a;|^ft^(r)cr. V*+(r) 
+¥_{r)a* ■ V^'_(r)| , 



where vp = ^ is the Fermi velocity (which is hence 
forth set to unity, along with h), ^j.{r) = (aj_(r), &^j_(r)) 
and Pauli matrices are defined as ti = {ax,<Jy) & cr* = 
(axj—cTy) Above Hamiltonian accurately describes 

the dynamics of electrons in an ideal flat graphene, free 
from any disorders or dislocations. The effects of disorder 
or deformation are captured by Hq, since deformations 
would lead to change in the hopping amplitudes locally, 
thereby giving rise to Sti (r) . It is shown in Ref . [1, [13] , 
that if Sti{r) is smooth over lattice spacing scales, then 
the two Dirac points are not coupled by disorder, and 
hence one can write : 

Hg^ y"d2^{*^+(r)a-a(r)*+(r) 

-^''Lir)(T* ■ a(r)^'_(r)| , 

where a = (ai,a2) is deformation induced gauge field, 
and is defined as ai(r) = ^ti(r) — ^{St2{r) + St^lr)) and 

02 (r) = ^{St2{r) — (5^3 (r)). It is very interesting to 
note, that the disorder couples to fermion field as a gen- 
uine U(l) gauge field, since SUs are defined on the links 
connecting the two sublattices, and unlike other cases 
where the effects of disorder can be captured by averag- 
ing method or replica method. The appearance of disor- 
der as an effective gauge field, is not astonishing. It is 
known that effects of dislocations (like presence of an odd 
pentagon/heptagon in the honeycomb lattice), local cur- 
vature and phonon-electron interactions can be treated 
as gauge fields d, [2lj . Interestingly, certain phonons in 
graphene couple to matter current as ip"f^A^^, rather 
than coupling as '00'/' j giving rise to some interesting 
effects that were seen experimentally by Raman spec- 
troscopy HH-ili. 

In analogy with Dirac equation for massless particles, 
one can unify the two spinors by working in a re- 
ducible representation of Clifford algebra, and defining 



(density) now reads: 



(3) 



where i,j run from 1,2 (sum over repeated indicies is 
implied here and hence forth). In above Hamiltonian, 
we have introduced explicitly coupling constant g, which 
specifies the strength of the fermion-gauge field interac- 
tion, and we have used the standard Dirac representation 
for 7 matrices as : 
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(z- 1,2,3), 
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(4) 
(5) 
(6) 



and (Tj's are the Pauli matrices. Taking into account, 
the non-dynamical component oq, to ensure that the 
ground state to start with, has zero charge density, the 
Lagrangian corresponding to above Hamiltonian, with a 
mass term, may be written in relativistic form as: 



C = — m - 



(7) 



which has an emergent Lorentz symmetry. Interest- 
ingly, the matter current, ipj^j^'y^ip coupled to the 
gauge field is conserved, even in the presence of a mass 
term. The corresponding conserved charge is Q ~ 

Jd^x\^{ala+ + blb+) - {ala^+blb^)Y 

In above Lagrangian, we have introduced a mass 
term, which does not follow from the hopping Hamilto- 
nian. The massless nature of the fermionic excitations in 
graphene in presence of a substrate have been confirmed 
conclusively from experiments @. However, recent com- 
putational studies by Drut & Lahde has attracted signifi- 
cant attention, who showed that free suspended graphene 
in vaccum, may well be an insulator with a gap, in a 
resonanting valence bond ground state [23 - |27| . Earlier, 
on the grounds of symmetry and renormalization group 
analysis, it was shown that at low energies, graphene may 
break chiral symmetry and can become gapped [28j . Un- 
fortunately, the experimental study on freely suspended 
graphene sheets have been difficult and hence the nature 
of its ground state is hence unclear. Also, it has been 
shown that, graphene can infact be made gapped, that 
too by several ways, for example by inducing local po- 
tentials [27|. Hence, consideration of gapped graphene 



spectrum is not an unphysical assumption. 

The external photon field, living in 3 -I- 1 D space time, 
interacting with fermions confined on a plane can not be 
correctly described by -JFf,^F^'•' {p., v = Q, 1,2), since 
it has only one degree of freedom, and not all compo- 
nents of physical photon field couple to matter. Follow- 
ing the arguments of Ref. [29j . we begin with photon 
field action / d'^x^Ff^yF^^'^ -\- j^A^^, where describes 
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the matter current confined to XY-plane. Integrating 
out z-coordinate in the above action, and using Greens 
function identity, we get the three dimensional action as 
: / (Px^F^,^F^''' ^j^A^', where F^, {^J^, = 0, 1, 2) 

describes 3D physical magnetic field Bz, in-plane com- 
ponents of electric field and Ey. Hence, the above 
Lagrangian suitably coupled to photon field reads: 



1 



(8) 



It needs to be pointed out here that, we have assumed 
velocity of light c equal to the Fermi velocity vp which is 
set to unity. This is actually not the case, but as shown 
in Ref. [29|, [3d| , this does not change the qualitative be- 
haviour of the theory. Also, notice that the coupling g 
is not dimensionless, as in QED^, but rather has dimen- 
sion of square root of mass, and the potential between 
two such charges is logarithmic, as expected. Next we 
integrate out the fermion fields, and find one loop con- 
tribution to effective action using derivative expansion of 
fermion determinant [sj, S3 , as : 



1 



-eff 
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^""^ 2-K\m\ 



As one can clearly see, the two gauge fields are cou- 
pled by a mixed Chern-Sinions term, which has a topo- 
logical nature [33]. Note that the pure Cherm Simons 
term for both gauge fields, in this case is absent, as 
rr(7^7'^7''7^7^) is the only nonvanishing trace of three 
gamma matrices in this representation, and gives rise 
to mixed Chern-Simons term. Also note that the sign 
of m appears as the coefficient of this mixed term, and 
hence absolute magnitude of the mass is irrelevant for this 
mechanism, and of course massless fermions will forbid 
this coupling. The Chern-Simons coupling, may lead one 
to conclude that parity must be violated in this interac- 
tion. However, by defining A± = A±a fields, the Chern- 
Simons term in above Lagrangian can be rewritten, mod- 
ulo a constant, as : e^^iypA^d'^ A'^ — e^^pAtd'^ A''_^ which 
has two topological terms with opposite signs; as a result 
it does not violate parity. It is straightforward to see, 
by integrating out a^, that the above coupling induces 

a parity conserving non topological mass M = '^^ j™^ 
for the physical electromagnetic field and leads to spon- 
taneous breaking of the Ue{^) symmetry of the elec- 
tric charge. The photon field becoming massive im- 
plies Meissner effect, a hallmark of superconductivity, 
where the static magnetic field exponentially dies down 
with distance from the boundary, with characteristic 
length scale, called penetration depth, which in our case 
is A = g-^^l^i . To make the comparison with usual 
BCS theory clear, we follow the technique given in Ref. 
[34], and define an auxiliary scalar field such that 



Hubbard- Stratonovich transformation, reads: 



C,,,^Un^\m\[d,<l>-£^A, 



(10) 



As we can clearly see, the above Lagrangian is in man- 
ifest London form, and all the phenomenological prop- 
erties of superconductivity would follow from here [35j ■ 
It is worth pointing out, that this is a non-BCS type of 
superconductivity, and is present only in the planar 2 -t- 1 
D world. 

Following the arguments given in Ref. [1^ [s^ , one can 
see that the correspondence between and (f> ceases to 
be one-to-one when vortex excitations are present in a^, 
and the above Lagrangian has to be corrected to consider 
this effect. It can be done by substituting dfj.4> by 9^0 — 
iip*dfj,ip, where ip*ip ~ 1. Defining vortex current as: 
jvort ^ {2TTi)-h^,^xd''{f*d^(p), allows us to write the 
Lagrangian : 



-evr^mse^.AJ™'-*-''^^''^^ - ^M!f^,lF^^ (11) 

Notice that the last term is topologically trivial, and 
is responsible for superconductivity. At low tempera- 
tures, the vortices are bound and at length scales larger 
than the binding scale, photon field sees vortex neu- 
tral current i.e., J™''* « 0. Hence, there is no contri- 
bution from the vortex terms in the above Lagrangian 
and one has a superconducting phase. The vortex 
bound states have finite energy, and at a finite temper- 
ature Tbkt the vortex unbinding takes places, known 
as Berezinskii-Kosterlitz-Thouless phase transition in the 
literature (stI . Issj . Considering a finite vortex current, 
one can show that the vortex contribution to photon 
field is •^l™^g p^^^piJ-t^ ^ which precisely cancels the con- 
tribution from non- vortex part, hence destroying super- 
conductivity. So, the superconducting-to-normal phase 
transition occurs at temperature T = Tbkt, which is 
given by following transcendental equation: 



tanh 



\m\l3_ 



BKT 



Pbkt — 



1 



67r31 



(12) 



^niypd'^a'', so that the above Lagrangian after a 



This equation, when numerically solved, gives Tbkt = 
■1^, where a — 0.1037312 [s^. As one can clearly see, 
Tbkt oc \m\. 

We have found that, disorder can couple to massive 
fermionic quasiparticles in graphene as a gauge field, 
to give rise to a special kind of superconductivity. As 
pointed out by Sasaki and Saito, such a gauge field can 
also rise from certain phonon modes (2l| . In usual pairing 
mechanisms, which describe low temperature supercon- 
ducters, an effective four fermi interaction leads to cre- 
ation of fermion-fermion pairs, which successively Bose 
condense to give a new superconducting ground state, 
which has lower symmetry than the Lagrangian. Hence, 
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in such case, superconducting state is a new vacuum, 
which is unitarily inequivalent to normal perturbative 
vacuurn, which has free fermionic excitations asymptot- 
ically [40'|. Further, the superconducting vacuum is a 
condensate of fermion pairs, and is a macroscopic object. 
Certain excitation of this condensate, so called Ampli- 
tude Mode, have been seen experimentally, establishing 
unambigiously its existence. In sharp contrast to this, the 
current proposal for superconductivity in graphene, is not 
described by any kind of pairing whatsoever. We show, 
that the perturbative ground state of gapped graphene, 
in presence of a gauge field which may arise from dis- 
order, exhibits superconductivity. This directly implies 
absence of any amplitude mode, since there is no conden- 
sate in this model. This feature makes this type of su- 
perconductivity experimentally distinguishable from the 
conventional ones. 

An interesting and very promising possibility of ob- 
serving this type of superconductivity is in ultra cold 
fermionic atoms trapped in two dimensional optical lat- 
tices having hexagonal honeycomb structure. It is re- 
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